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ABSTRACT 

The gravitational lens 0957+561 was monitored with the Very Large Array 
Q\ \ from 1979 to 1997. The 6 cm light curve data from 1995-1997 and the 4 cm 

data from 1990-1997 are reported here. At 4 cm, the intrinsic source variations 
occur earlier and are twice as large as the corresponding variations at 6 cm. 
The VLBI core and jet components have different magnification factors, leading 
7— i \ to different flux ratios for the varying and non-varying portions of the VLA 

light curves. Using both the PRHQ and Dispersion statistical techniques, we 
determined the time delay, core flux ratio, and excess non-varying B image 
| flux density. The fits were performed for the 4 cm and 6 cm light curves, both 

individually and jointly, and we used Gaussian Monte Carlo data to estimate 
68% statistical confidence levels. The delay estimates from each individual 
wavelength were inconsistent given the formal uncertainties, suggesting that 
there are unmodeled systematic errors in the analysis. We roughly estimate the 
systematic uncertainty in the joint result from the difference between the 6 cm 
and 4 cm results, giving 409 ± 30 days for the PRHQ statistic and 397 ± 20 days 
for the Dispersion statistic. These results are consistent with the current optical 
time delay of 417 ± 3 days, reconciling the long-standing difference between the 
optical and radio light curves and between different statistical analyses. The 
unmodeled systematic effects may also corrupt light curves for other lenses, 
and we caution that multiple events at multiple wavelengths may be necessary 
to determine an accurate delay in any lens system. Now that consensus has 
been reached regarding the time delay in the 0957+561 system, the most 
pressing issue remaining for determining H is a full understanding of the mass 
distribution in the lens. 
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(0957+561) 



1 Haverford College, Haverford, PA 19041, dhaarsma@haverford.edu 

2 Department of Physics, 37-607, Massachusetts Institute of Technology, Cambridge, MA 02139 
3 Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 



-2- 



1. Introduction 

The time delay between multiple gravitationally lensed images can be used to measure 
the distance of high redshift objects, and thus is a useful estimator of the Hubble parameter, 
Hq. After many years of monitoring the lens 0957+561, the time delay estimates for this 
system are finally converging on an accepted value. Groups monitoring 0957+561 at optical 
wavelengths have detected a sharp variation in each image, and have found the optical 
delay to be 417 + 3 days flKundic et al. 1995 , 1997; Uscoz et al. 1997 ; Schild fc Thomson 



|1997|) . Given the long controversy over the value of the delay (for a history, see Table 1 of 
Haarsma et al. 1997, hereafter Paper 1), it is important that the optical measurement be 
confirmed at radio wavelengths. The MIT radio astronomy group has monitored the source 
at radio wavelengths from 1979 to 1997 and the final light curve data and time delay results 
are reported here. 



2. Observations 

Observations have occurred monthly at the National Radio Astronomy Observatory 
(NRAO) Very Large Array radio telescope (VLA)Q since 1979 at 6 cm and since 1990 at 
4 cm. The monitoring ended in December 1997. All of the data were reduced in the manner 
described in Paper 1 and Lehar et al. (1992). To determine the flux densities of the point 
images, it was necessary to subtract the extended structure in the field. At both 6 and 
4 cm, this subtraction was difficult in the most compact VLA array, D, thus there are gaps 
in the light curves for 4 months of every 16 month cycle. In addition, some observations in 
other VLA arrays were excluded due to bad weather or poor subtraction of the extended 
structure. When the observations were made in a combination or non-standard array 
configuration, the data were analyzed according to the next largest standard configuration 
(A, B, or C). 

The 6 cm data through December 1994 were presented in Paper 1. The 
remaining 6 cm data and all of the 4 cm data are given in Tables 1 and 2, and 
plotted in Figure 1. The light curve data are also available electronically through 
[http : / /space . mit . edu/RADIO/papers . html| . There are a total of 147 points in the 6 cm 



light curve, and 58 points in the 4 cm light curve. At 6 cm, the flux density of the B image 
increased in 1995, following the A image increase in 1994. The current 6 cm feature has 



4 The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under 
cooperative agreement with the National Science Foundation. 



- 3- 



lasted longer than the similar feature around 1989-1991, but the A image is now declining. 
At 4 cm, the quasar is twice as variable as at 6 cm (as a percentage of average flux density). 
Also, the variations in the 4 cm light curves occur earlier than the corresponding features 
at 6 cm. Both of these characteristics are consistent with multi-wavelength models and 
other observations of AGN variability (e.g. [Marscher fc Gear 19851 ; IStevens et al. 199^ ) . 
The well-sampled increase and decrease at 4 cm in 1994-97 has helped significantly in 
determining the radio time delay. 



3. Free parameters in the Light Curves 

When fitting for the time delay between the images, the difference in magnification 
between them must be properly taken into account. In past analyses of lensed light 
curves (including Paper 1), only two parameters were used in the fit: the time delay and 
a single flux ratio. Conner, Lehar, & Burke (1992), however, have pointed out that the 
magnification varies rapidly along the B image, causing the VLBI core and jet components 
to have different flux ratios, with the core ratio being larger. At the resolution of the VLA, 
the beam includes both the core and the jet, and thus the flux ratio of the VLA light 
curves -Rvla is a composite of the core and jet values. The VLA light curve is the sum of 
the jet (which is constant in time), and the core (which has both variable and constant 
components). There are then four physical parameters: the time delay r, the flux ratio 
of the core -R CO re = B core /A COTe , the flux ratio of the jet i?j Ct = Bj et /Aj et , and the amount 
of flux density due to the jet vs. the core, i.e. B(t) = B core (t) + B- iCt for the B image. 
Note that the core can contain both constant (DC) and variable (AC) components, i.e. 
S core = -B cor e,DC + Bac- Also, the DC part of the light curve is due to both the core and the 

jet, i.e. B DC = -Bcore.DC + -Bjet- 

Press & Rybicki (1998) discuss these issues in the context of the optical light curves 
of 0957+561. They point out that the amount of constant flux due to the core (-B C ore,Dc) 
is impossible to determine, since we may have not yet seen the variable part of the core 
(B\c) go to zero. They show that on a fundamental level there are only three measureable 
parameters in a pair of lensed light curves, which may be cast as: the time delay r, the 
core flux ratio R COTe , and the extra constant flux in the B image that does not occur in the 
A image, 

c=p A DC (1) 

-Kac 

( Press fc Rybicki 1998|) , where Rac is the flux ratio of the variable component. It is useful 
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to write c in terms of the core and jet components of the radio images as 

~ (-Ajet + A:ore,Dc) , (2) 



-Bjet + B 

core, DC 



-Rcore 

and therefore 

= Sjet f - , (3) 

\ -fl-core -^jet / 

where the DC core components cancel out. The value of c can thus be estimated from the 
values of £?j et , -R CO re, and Rj et . Since in the case of 0957+561 we have Rj et < R CO re flConner 



Lehar, fc Burke 1992| ), the value of c must be negative; i.e. the A curve has a larger amount 



of constant flux than the core-ratio corrected B curve. 

The values of several of the above parameters can be estimated from observations 
without doing time delay fitting. Garrett et al. (1994) compiled the information on the 
core flux ratio from VLBI and optical observations, and found the weighted average of these 
estimates to be -R cor e = 0.75 ± 0.02. Also, the faintest portions of the VLA light curves 
set upper limits on the jet flux density, i.e. £?j et ^ 21 mJy at 6 cm, and Bj et <^ 15 mJy 
at 4 cm. A better estimate of B ]Ct can be obtained by comparing coincident VLBI and 
VLA observations. The VLBI observations give the core flux density at a particular 
epoch, which can be subtracted from the VLA flux density to obtain the VLA jet flux 
density. Campbell et al. (1995) report VLBI observations at 6 cm on 1987 Sep 28 and 
1989 Sep 26, and by comparing these to VLA observations occurring on the same days we 
find -Bj et = 11.1 ± 0.4 mJy and Rj et = 0.63 ± 0.03. The values for -R cor e, -Rjet, and -B CO re can 
be combined using equation [| to find cq = —2.7 ± 0.8 mJy. 

The above estimates are all for the 6 cm light curves. At 4 cm there are no coincident 
VLBI/ VLA observations, so we can not make similar estimates. The value of c is different 
at 6 cm and 4 cm due to the difference in Bj et ; note that the ratios R iet and -R CO re are the 
same for the two bands. For a synchrotron spectrum, 5j et will be smaller at 6 cm than 
4 cm, and thus we expect \c^\ to be smaller than \cq\. 



4. Time Delay Analysis Methods 

To fit for the three parameters r, -R core , and c (described in we used the PRHQ 
statistic (Press, Rybicki, & Hewitt 1992a, 1992b; |Rybicki fc Press 1992| ; incorporating the 
modifications of |Rybicki fc Kleyna 1994] ; Press fc Rybicki 1998j ), and the Dispersion statistic 
( Pelt et al. 1994j , 1996), which were described in Paper 1. We used linear units (mJy) rather 
than the logarithmic units defined in Paper 1. The discrete correlation function ( [Lehar et al. 



1992| ) did not find a strong correlation in the 4 cm light curves, so that statistic was not 
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used here. Gaussian Monte Carlo data were made as described in Paper 1, but now with 
the four physical parameters r, -R cor e, Rjet, and Bj et . Five hundred Gaussian Monte Carlo 
data sets were used to estimate the 68% confidence intervals on the results for the real light 
curves, where the fitted c values were compared to the input parameters using equation 
The pseudo-jackknife test from Paper 1 was used to test the stability of the result to the 
removal of individual points. 

To determine whether neglecting the difference between the core and jet flux ratio 
caused an error in our previous analysis, we applied the two dimensional fit (for r and 
-Rvla, as in Paper 1) to the Gaussian Monte Carlo data made with four parameters. The 
resulting fitted-minus-true values did not show a significant bias (to long or short delays, 
for example), but did show an increase in scatter about the true delay. We found that the 
error in the delay increased monotonically with Bj et , from roughly 20 days for B ]et = mJy 
to roughly 100 days for B ]et ~ 11 mJy (the value for the real light curves) when using 
the PRHQ statistic and the 6 cm Monte Carlo data. The same test with the 4 cm Monte 
Carlo data, and with the Dispersion statistic at both 4 cm and 6 cm, revealed a similar 
but somewhat milder effect, with the delay error at least doubling between small and large 
values of B^ et . This dependence on the amount of flux density in the jet component may be 
one cause of the inconsistency in delay estimates over the years, and we caution that fitting 
for only two parameters may introduce significant errors. 

In addition to analyzing the two wavelengths individually, we also fitted for the 
parameters using both wavelengths at once. The Dispersion and PRHQ statistics are 
both easily modified for this by minimizing the sum of the statistics from each wavelength 
(see Press et al. 1992b), and fitting for the parameters r, -R core , c 6 , and c 4 . Monte Carlo 
analysis was also done for the joint data, where the 6 cm and 4 cm Monte Carlo sets were 
constructed with the same set of [t, R core , Rj e t, -Bjetj- 

The covariance model (Paper 1; Press et al. 1992a) used for the PRHQ statistic was 
found by an iterative procedure on the individual light curves. First, measurement errors 
of 2% were assumed and used to make point estimates for the structure function, and 
then fitted to an exponential in the lag range of 100 to 700 days. This structure function 
was then used to determine the PRHx 2 value for the light curve, and the measurement 
errors were adjusted until PRHx 2 equaled the degrees of freedom. Then the process was 
repeated for the new measurement error value. Iterations stopped when the square root of 
(PRHx 2 /degrees of freedom) changed by less than 1% when PRHx 2 was calculated with 
the measurement error of the previous iteration. At 6 cm, the covariance model found was 

V{T) = 1.673 x 10~ 4 T L606 mJy 2 , (4) 

with measurement errors ca = 1.82% and e# = 2.34%, where T is the time lag between two 
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points on the curve. At 4 cm, the fitted covariance model was 

V{T) = 3.174 x io- 4 T L633 mJy 2 , (5) 

with measurement errors ca = 1.67% and e# = 2.19%. 

The sharp feature in the B image at 6 cm in Spring 1990 is statistically inconsistent 
with the rest of the light curve (see Paper 1), and, given the short delay found from the 
optical and 4 cm light curves, the A image shows that the feature is not intrinsic to the 
source. Thus we expect that the results with the points removed will be more accurate than 
the results for the full 6 cm curves. The 6 cm analysis was done both with and without 
the four points (1990 March 15, April 10, May 7, and May 23); the light curve without the 
points will be denoted 6*cm. 

5. Results 

The main results of the time delay analysis are shown in Tables 3 and 4. These and 
other aspects of the results are worth discussion in this section. First, the Monte Carlo 
analysis showed that the use of the three parameter fit (r, -R cor e, c) made the time delay 
confidence interval independent of B seXi (rather than increase with Bj et as happened in the 
two parameter fit described in §f|). Next, the pseudo-jackknife test (see Paper 1), using 
either the PRHQ or Dispersion statistic, showed that removal of an individual point from 
the light curve general caused a change in the delay that was much smaller than the 
confidence interval, with a few important points in the curve (typically during rises or falls) 
causing a change at about the amount of the confidence interval. Thus, the delay estimate 
in all cases is stable under the removal of individual points. Note also that the removal of 
the four Spring 1990 points (6*cm vs. 6 cm in Tables 3 and 4) never caused a change in the 
fitted parameters of more than the confidence interval. 

The relationship between the values of c and -R cor e for a given delay is worth pointing 
out. Figure 2 shows the PRH% 2 statistic as a function of c and -R cor e for the 6*cm light 
curves, with the delay fixed at the best fit value. The surface is a diagonal trough, with 
the location of the minimum poorly constrained along the bottom of the trough. The 
trend is such that a larger value of -R cor e requires a more negative value of c, as expected 
from equation |3|. The 4 cm light curve has a similar PRHx 2 surface, and the Dispersion 
surface shows a similar but even more pronounced effect. Thus, if either -R cor e or c is poorly 
constrained by the light curves, the other parameter will also be poorly determined. 

There is a significant bias in c and R COT e f° r the Dispersion statistic, which is related to 
the interdependence of these parameters. The fitted-minus-true values from the Dispersion 
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analysis of the Gaussian Monte Carlo data showed a pronounced asymmetry and bias in 
both c and -R cor e (but was nearly symmetric and un-biased in delay). Tests showed that the 
bias in c and -R cor e was somewhat reduced for S ~ 25 days (where 5 is a weighting parameter 
in the Dispersion statistic, see Paper 1). We continued to use 6 = 60 days (as done by 
pelt et al. 1996|) , since that value had the narrowest distribution for the fitted-minus-true 
delays. Thus, only the time delay results are given in Table 4. 

The PRHQ statistic produced a symmetric distribution in c and -R CO re, and the results 
are listed in Table 3. Note that the fitted value of c is more negative when -R cor e is larger, 
as mentioned above. The analysis of the individual wavelengths found somewhat different 
values of -R cor e, and that c 4 is more negative than c 6 contrary to the predictions of §|3|. 
The joint analysis of the two wavelengths, however, used all of the available information 
to constrain the fit, and the resulting i? CO re, Cq, and C4 are in good agreement with the 
predictions of 

Turning now to the results for the delay, it is of interest that all of the delay estimates 
in Tables 3 and 4 are much smaller than the value of 540 days found by Press et al. (1992b) 
using the first 80 points in the 6 cm light curve. As explained in Paper 1, this change in the 
delay estimate is due entirely to the addition of new features to the light curve. If the first 
80 points in the curve (with or without the four Spring 1990 points) are fitted for the three 
parameters, the Dispersion statistic finds a delay of roughly 550 ± 35 days, and the PRHQ 
statistic (using the above covariance model, eq. |4j) finds a delay of roughly 525 ± 25 days 
(the confidence intervals were estimated from Monte Carlo analysis of 100 light curves). 
Therefore the change in the delay estimate between the first 80 points and the current 147 
points occurs for both statistical methods, for fits with two or three parameters, and for 
a variety of covariance models. It does not, however, occur in our Gaussian Monte Carlo 
data. Comparing delay estimates from the first 80 points and the full curves in the Monte 
Carlo data, we find that 99% of the sets have a difference in delay less than the 75 day 
difference seen in the real light curves. Thus the Monte Carlo light curves must still be 
missing some characteristic of the real data. 

A related issue is the difference in the delay estimates for the two wavelengths (for 
PRHQ, r 6 = 452i}| days vs. T4 = 397 ± 12 days, a difference of about three confidence 
intervals). The time delay between lensed images should be completely independent of 
wavelength, and indeed the optical and radio estimates come remarkably close. This effect 
also does not occur in the Gaussian Monte Carlo data, which were constructed such that 
the 6 and 4 cm data have the same set of [r, R core , Rjet, -Bjet]- When applying the PRHQ 
statistic, we found that 99% of the Monte Carlo curves had a difference of (\tq — T4J) that 
was smaller than the 55 day difference in the real data; similarly, the Dispersion statistic 
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found that about 80% of the data sets had a smaller delay difference. 

Since both of these effects (the significant change in the delay estimate as features are 
added to the curves, and the significant difference between the two radio wavelengths) are 
not seen in the Monte Carlo data, there must still be some systematic effect that has not 
been taken into account in the creation of the Monte Carlo data or in our analysis. One 
source of the systematic error may be interstellar scintillation, which can cause variability 
at a level not much larger than our observational error of 2%, perhaps creating small 
features such as the discrepancy between the 6 cm A and B images in early 1985 (see 
Figure 4). Although individual features of a few percent are difficult to identify, they may 
cause a significant bias in the delay estimate if they occur at crucial times in the light 
curves. We note that microlensing could cause similar low level systematic effects in the 
optical light curves QSchild fc Smith 199 1| ; [Schmidt fc Wambsganss 1998| ). Since it is beyond 



the scope of this paper to model these systematic effects, we make a rough estimate from 
the difference between the 6 and 4 cm delay estimates, giving a systematic uncertainty of 
roughly ±30 days for PRHQ and ±20 days for the Dispersion. Note that this uncertainty 
is still less than 10%, and is only one factor contributing to the error in H (see §6). Note 
also that the optical time delay estimate has a smaller error (only l%, [Kundic et al. 1997Q 



primarily because the source varies much more rapidly at optical wavelengths than radio 
wavelengths. We caution others monitoring gravitational lenses that in order to determine 
an accurate time delay it is preferable to have light curves with multiple features at multiple 
wavelengths, since the estimate of the time delay based on a single feature at a single 
wavelength could easily be corrupted by these low level systematic effects (as the 0957±561 
6 cm curves were after the first 80 observations). The 0957 radio light curves will continue 
to be a useful data set for studying systematic effects and time delay analysis techniques. 

Figure 3 shows the PRHQ and Dispersion statistics as a function of delay for the joint 
analysis of the 4 cm and 6*cm curves, with the values of c 6 , c 4 , and R CO re fixed at the best 
fit values. Figure 4 shows the aligned light curves at the two wavelengths with the PRH 
optimal reconstruction (see Press et al. 1992a and Paper 1). 



6. Conclusions 

Since our last report (Paper 1), the B image has increased at 6 cm, and the A image 
has entered a slow decline. The 4 cm curves, given here for the first time, are highly variable 
and give additional features to aid in determining the time delay. To take into account the 
difference in magnification of the core and jet components of each image, we have fit for 
the delay, core flux ratio, and the excess flux density in the B image (as defined by Press 
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& Rybicki 1998). The delay estimates found from the wavelengths individually disagree 
by a few confidence intervals, indicating that there are systematic effects not modeled in 
our analysis. The delay estimates found from the joint analyses of both wavelengths were 
409 ± 30 days for PRHQ and 395 ± 20 days for the Dispersion, where the uncertainty is 
based on a rough estimate of the systematic error. Both results are consistent with the 
delay estimated from optical monitoring (417 ± 3 days, Kundic et al. 1997|) , and thus we 



now have good agreement for the value of the delay from both statistics in both the radio 
and optical light curves. Consensus has finally been reached on the value of the delay for 
gravitational lens 0957+561. 

This measurement of the delay can now be used to answer cosmological questions. 
The Hubble parameter, however, depends not only on the delay but also on the lens model 
and the galaxy velocity dispersion. Using the SPLS model of Grogin & Narayan (1996a, 
1996b), the recent Keck velocity dispersion measurement of 279 km s -1 flFalco et al. 1997| ). 
and a time delay of 409 days, we obtain ifo=67kms _1 Mpc~ 1 . In fitting this model, Grogin 
& Narayan used the ground-based optical position of the lensing galaxy as a constraint, 
rather than the more precise VLBI position. Since then, the HST optical position of 
the lensing galaxy has been found to agree with the VLBI position ( [Bernstein et al. 



1997| ). The modelers point out flGrogin fc Narayan 1996 b| ) that if the VLBI position is 
used as the model constraint, their model fit is very similar to that of Falco, Gorenstein, 
& Shapiro (1991); for the same delay and velocity dispersion, the Falco et al. model gives 
H =A1 kms _1 Mpc _1 . Thus, the change in the position of the lensing galaxy causes a change 
in the estimate of Hq well beyond the statistical error, apparently in contradiction to the 
conclusions of Kundic et al. (1997) regarding the robustness of the Hq determination. New 
modeling work must be done which incorporates the improved galaxy position, the recent 
observations of the cluster mass distribution ( Fischer et al. 1997| ), a careful treatment of 
the systematic errors in the velocity dispersion ( [Romanowsky fc Kochanek 1998| ), and the 



recently reported structure at x-ray ( |Chartas et al. 19"98| ), optical ( [Bernstein et al. 1997| ) 



and radio ( |Avruch et al. 199?] [Harvanek et al. 1997| ) wavelengths. These new observations 



will allow an improved fit of the model to the data and provide a more accurate measure of 
the Hubble parameter. 
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Fig. 1. — The complete 6 cm and 4 cm light curves of gravitational lens 0957+561. The A 
image data are shown as triangles and the B image as circles. The 4 cm A image has been 
shifted up by 8% to avoid overlap with the B image. 
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Fig. 2. — PRHQ for the 6*cm light curves, as a function of c and R corc . The delay is fixed 
at 452 days. The minimum is Q = 125.6 at c = —1.47 and R CO re = 0.731. Contours start at 
Q = 130 and increase by 10 to Q — 380. 
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Fig. 3. — Joint time delay analysis of the 6* cm and 4 cm light curves. The PRHQ statistic 
is shown as a solid line with a minimum at 409 days. The Dispersion statistic is shown as 
a dashed line with a minimum at 395 days. The values of -R CO rc, c&, and C4 were set at the 
best fit values for purposes of plotting the statistic vs. delay. The vertical axes were scaled 
such that the PRHQ and Dispersion confidence intervals were the same height. 
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Fig. 4. — The 6 cm and 4 cm light curves combined at at r = 409 days, -R CO rc = 0.753, 
c 6 = —2.35, and C4 = —2.09, shifted to the time and flux density of the A image. The A 
image data are shown as open circles and the B image as solid circles. The one sigma width 
of the PRH optimal reconstruction is shown as a pair of lines. 
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Table 1. 6 cm Light Curve Data 
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9886.57 


D^A 


OCT 1 O 

35.13 


22. o9 


1995 Jun 23 


9892.50 


A 


OP oo 

36.82 


OO P'V 

26.01 


1995 Jun 28 


9896.53 


A 


OT ( \ O 

37.02 


O O <7C 

22.75 


1995 Jul 08 


9907.23 


A 


or r\fy 

35.97 


O C OO 

lb. 61 


1995 Jul 21 


9919.51 


A 


36.13 


a no 
24. 92 


1995 Aug 07 


9937.34 


A 


op no 

36.02 


O A PC\ 


1995 Sep 01 


9962.33 


A 


O r 1/1 

35.14 


a Tr 

24.75 


1995 Sep 09 


9970.14 


A^AnB 


37.12 


O A C\A 

24. U4 


1995 Sep 15 


9976.10 


AnB 


OO 91 

08.0I 


O A 

24.oo 


1995 Sep 23 


9984.17 


AnB 


r. nn 

00.99 


O A AQ 

24. 4o 


1995 Sep 30 


9991.17 


AnB 


OP OP 

36.36 


25.90 


1995 Oct 10 


10001.17 


B 


OP PO 

36.68 


OP 1 O 

26.13 


1995 Oct 27 


10018.20 


B 


OP r n 

36.50 


or op 

25.36 


1995 Nov 09 


10031.06 


B 


or on 

35.39 


A no 

24.08 


1995 Dec 26 


10077.95 


B 


35.66 


O A 

24.65 


1996 Jan 26 


10108.83 


BnC 


00 nn 

33.99 


or 10 

25.13 


1996 Feb 05 


10118.77 


BnC 


or 10 

35.18 


O /I r ^7 

24.57 


1996 Feb 26 


10139.67 


C 


O C AO 


a r\T 
24. (J7 


1996 Mar 04 


10146.68 


C 


OP E E 

00.00 


24. lo 


1996 Apr 05 


10178.56 


c 


OA Ol 

34.31 


r nr 

25.06 


1996 Apr 25 


10198.66 


c 


rr r\r\ 

35.09 


or AT 

25.07 


1996 Jun 11 


10246.48 


CnD 


oa nn 

34.00 


rir -1 A 


1996 Oct 19 


10376.05 


A 


35.30 


25.20 


1996 Nov 10 


10397.95 


A 


34.62 


25.36 


1996 Dec 26 


10443.98 


A 


35.17 


24.74 


1997 Jan 10 


10458.89 


A 


34.81 


24.86 


1997 Feb 26 


10505.82 


B 


33.62 


24.38 


1997 Mar 19 


10526.69 


B 


33.95 


25.53 


1997 Apr 10 


10548.65 


B 


33.92 


25.90 


1997 May 11 


10579.65 


B 


34.20 


25.06 


1997 Jun 22 


10622.38 


BnC 


32.38 


24.60 


1997 Jul 11 


10641.43 


C 


33.28 


24.49 


1997 Sep 22 


10714.17 


C 


33.00 


23.73 


1997 Oct 06 


10728.18 


CnD 


33.06 


22.32 
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Table 1 — Continued 



Calendar Day a Array Flux Density (mJy) 



Date A Image B Image 



a Julian Day - 2,440,000.0 
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Table 2. 4 cm Light Curve Data 


Calendar 


Day a Array Flux Density (mJy) 


Date 


A Image B Image 



1 QQD 

X C/C/U 


Ort 04 


81 fiQ 
o luy. 


99 

. ZfZ; 






9S 

ZjO . 


87 


91 

_ X . 


Qfi 

. C/U 


1 QQD 


Nov Q1 


81 Q7 

O X U \ , 


06 

UU 




c 


23 

ZjO . 


71 

. 1 X 


99 


57 

■ O I 


1 QQD 

x u u u 


Der 1 3 


8238 


8Q 




c 


22 

Z Z . 


7 s ) 


99 

Z Z . 


4Q 


1 QQ1 

x uu X 


Tni m 


8448 

(J tit: (J , 


40 

T:U 




A 

A 


22 


28 


1 7 

X l . 


Q8 


1 QQ9 

X C/C/Zj 


Tan flfi 

J dll uu 


8fi97 


Q7 




R 


ZjO . 


■ OO 


xo . 


1 9 

. X Z 


1 QQ? 

X C/C/Zj 


Fph 04 

1 CU Ut: 


OUuU. 


80 




Rnf! 


99 

ZfZ; . 


Q^ 


1 ^ 

xo . 


^4 


1 QQ2 

X C/C/Zj 


Fph 9Q 


8fi81 

OUO X , 


74 

■ I t: 




f! 


94 

Zjt: . 


08 

.UO 


1 ^ 

xo . 


79 


1 QQ2 

X U U Zj 


Mar D7 

1V10j± U 1 


8688 


67 

. U I 




C 


93 

ZjO . 


Q6 


1 ^ 

X o . 


. OO 


1 QQ2 

X U U Zj 


Anr 1 R 


8730 


60 

. UU 




c 


94 

ZjT: . 


30 

. OU 


1 ^ 

X o . 


8 s ) 

. OO 


1 QQ2 

X C/C/Zj 


IVXclj uo 


O 1 TlU . 


fiO 




f! 

V 


94 

ZjT: . 


Q8 


1 ^ 

xo . 


88 

-OO 


1 QQ2 


Nov 1 1 

1 N \J V X J- 


8Q38 


0Q 




A 


95 

ZjO . 


^7 

-O I 


1 ^ 

xo . 


61 

. U X 


1 QQ2 


Dpr ID 


8Q66 

O <-/UU . 


Q7 




A 


95 


1 Q 

. X u 


1 ^ 

X o 


34 

. Ot: 


1 QQ3 

X <_/ <_/ u 


Fph 0^. 

X CU UU 


Q023 


78 

. 1 O 




AnR 


95 


99 

. _ Z 


1 ^ 

xo . 


68 

. UO 


1 QQ3 


Mar 21 

1VJ.C11 Z J- 


Q067 


64 

. Ut: 




R 


95 


^8 


16 

X u . 


91 


1 QQ3 

X U u o 


Anr DQ 

AUI U<_/ 


Q086 


67 

. U 1 




R 


9^ 

ZjO . 


7 s ) 

. 1 O 


16 

X u . 


78 

. 1 O 


1 QQ^ 

X C/C/O 


Mav 1 R 
ivxcty xo 


Q19fi 

C/ XZjU- 


48 


R 


—►Rnf! 


2^ 

ZjO . 




1 Pi 

xu . 


^4 


1 QQS 

X C/C/O 


Tnl 9 c i 

U U.1 ZjO 


Q1 Q4 

C/ X C/t:. 


21 

. Z; X 




f! 


27 


89 


1 7 

X I 


1 ^ 

. X O 


1 QQ3 

X U u o 


Aiip- 96 

A Lift ZjU 


Q996 

<_/ZjZjU . 


26 

. ZAJ 




C 


28 

ZjO . 


09 

. UZj 


1 6 

X u . 


6 s ) 

. UO 


1 QQ4 


Mar D4 

1VJ.C11 Ut: 


Q41 F5 


73 

1 U 




A 


31 

O X 


34 

. Ot: 


1 7 

X i 


71 

. f X 


1 QQ4 

X C/C/t: 


Anr 1 1 

iijJI X X 


Q4^S 

C/t:00. 


fi8 

. UO 




A 

A 


S1 

O X . 


1 4 


1 7 

X I . 


1 7 

. X l 


1 QQ4 

X C/C/t: 


Mm/ 07 
ivxdv u i 


Q47Q 

C/t: 1 C/. 


. UO 


A 


— > AnR 


O X 


S1 

.O X 


1 7 

X 1 


80 

.OU 


1994 


Jun 25 


9528. 


.52 




B 


31 


.15 


18. 


.58 


1994 


Jul 06 


9540. 


.42 




B 


30 


.83 


18. 


.56 


1994 


Aug 18 


9583. 


.28 




B 


31 


.78 


18. 


.62 


1994 


Sep 08 


9604. 


.27 




B 


31 


.65 


19. 


.42 


1994 


Oct 10 


9636. 


.18 




BnC 


31 


.04 


18. 


,97 


1994 


Nov 07 


9664. 


.08 




C 


31 


.75 


20 


,29 


1994 


Dec 08 


9694, 


92 




C 


31 


,32 


20 


.51 


1995 


Jun 23 


9892. 


.50 




A 


30 


.19 


20 


,99 


1995 


Jul 08 


9907. 


.23 




A 


30 


,49 


22 


.18 


1995 


Jul 21 


9919. 


.51 




A 


30 


,05 


21. 


.80 


1995 


Aug 07 


9937. 


.34 




A 


30 


.48 


21. 


.72 


1995 


Sep 01 


9962. 


33 




A 


29 


.94 


22 


.05 


1995 


Sep 09 


9970. 


.14 


A 


^AnB 


29. 


.91 


21. 


.02 
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Table 2 — Continued 



Calendar 


Day a 


Array 


Flux Density (mJy) 


Date 










A T 

A Image 


B Image 


1995 


Spn 1 i~i 

kJ \j VJ lv 


9976 


10 

. xw 


AnB 


30 76 


21 78 

til. 1 U 


1995 


Sen 23 


9984 


17 

■ X 1 


AnB 


29 45 


21 73 

Zj X . I KJ 


1 99 s ) 

X <_/ <_/W 


Sen 30 


9991 


1 7 

. X 1 


AnB 


29 19 

Zj . X <_/ 


22 ^4 

ZjZj . kj t: 


1995 


Ort ID 

V/l I X W 


10001 

X WWW X , 


17 

. X 1 


B 


29 41 

Zj u . T: X 


22.88 


1995 

X t_/ U W 


Oct 27 


10018 


20 

. Zj W 


B 


29 95 

Zj U . U KJ 


22.51 


1995 

X U U W 


Nov 09 


10031 

X W WO X , 


06 

. WW 


B 


28 55 

ZjU . UU 


21.81 


1995 


Dec 26 


10077. 


.95 


B 


28.34 


21.56 


1996 


Tan 26 

(J OjI J. — W 


10108. 


.83 


BnC 


27.02 


22.52 


1996 


Feb 95 


101 18 

X W X X o , 


77 


BnC 


27 19 

Zj I . X U 


21 96 

Zj X . U W 


1996 


Feb 26 

x \z v.f Zj w 


10139 

X W -L KJ 


67 

W 1 


C 


27 95 

Zj 1 . UKJ 


23 03 

ZjO .WO 


1996 


Mar 04 


10146. 


.68 


C 


28.11 


22.79 


1996 


Apr 05 


10178. 


.56 


c 


27.52 


22.89 


1996 


Apr 25 


10198. 


66 


Kj 


27.84 


23.15 


1996 


Oct 19 


10376. 


.05 


A 


25.49 


20.67 


1996 


Nov 10 


10397. 


.95 


A 


24.19 


20.64 


1996 


Dec 26 


10443. 


.98 


A 


24.73 


20.05 


1997 


Jan 10 


10458. 


.89 


A 


25.32 


20.13 


1997 


Feb 26 


10505. 


.82 


B 


24.53 


19.44 


1997 


Mar 19 


10526. 


.69 


B 


24.79 


19.56 


1997 


Apr 10 


10548. 


.65 


B 


25.23 


20.10 


1997 


May 11 


10579. 


.65 


B 


24.95 


18.65 


1997 


Jun 22 


10622. 


.38 


BnC 


24.71 


18.32 


1997 


Jul 11 


10641. 


.43 


C 


26.58 


18.31 


1997 


Sep 22 


10714. 


.17 


C 


26.97 


18.02 



a Julian Day - 2,440,000.0 
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Table 3. Results from PRHQ Statistic 



Light Curve 


l \ r 

Degrees of 
Freedom 


Q 


ni)ii 9 


Time Delay 

(days) 


Lore 
Flux Ratio 


c (mJy) 








6 cm 4 cm 






6 cm 


4 cm 


6 cm 


291 


165.5 


344 


459±H 


0.720l8:8S 






6*cm 


283 


125.6 


297 




n 7 oi +0.014 
u - ' ol -0.014 


_1 47 +o.56 




4 cm 


113 


15.7 


111 


397i^ 


7fi2 +0 - 026 

U./DZ_q Q31 




-2 44+ - 94 

z - 44 -0.77 


6 cm & 4 cm 




194.1 


352 115 


416±| 


u - '^-0.011 


z - uz -0.49 


i 7 o+0.37 
_i -'°-0.32 


6*cm & 4 cm 




148.6 


302 113 


409t^ 


A 7 CQ+0.011 


9 qr+0.45 


—2 nq +0 - 36 

z - uy -0.33 



*6 cm light curve with four points removed, see §4 



Table 4. Results from Dispersion Statistic 



Light Curve Dispersion Time Delay (days) 



6 cm 


0.608 


43011 


6*cm 


0.489 


41611 


4 cm 


0.259 


3831^ 


6 cm & 4 cm 


0.915 


3971H 


6*cm & 4 cm 


0.764 


3951H 



*6 cm light curve with four points removed, see 

§4 



